DO TH 99/22 
December 1999 



Massive Neutrinos in Astrophysics 

Lectures at the fourth national summer school 
"Grundlagen und neue Methoden der theoretischen Physik" 
31 August — 11 September 1998, Saalburg, Germany 

Georg G. Raffelt 

Max-Planck-Institut fur Physik 
(Werner-Heisenberg-Institut) 
Fdhringer Ring 6, 80805 Miinchen, Germany 

Email: raffelt@mppmu.mpg.de 

Lecture notes written by 

Werner Rodejohann 

Lehrstuhl Theoretische Physik III 
Universitat Dortmund 
Otto-Hahn Str. 4, 44221 Dortmund, Germany 

Email: rodejoha@dilbert.physik.uni-dortmund.de 



ABSTRACT 

An introduction to various topics in neutrino astrophysics is given for students with little 
prior exposure to this field. We explain neutrino production and propagation in stars, 
neutrino oscillations, and experimental searches for this effect. We also touch upon the 
cosmological role of neutrinos. A number of exercises is also included. 
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Preface 



Neutrino astrophysics is a prime example for the modern connection between astrophysics 
and particle physics which is often referred to as "astroparticle physics" or also "particle 
astrophysics." The intrinsic properties of neutrinos, especially the question of their mass, 
is one of the unsolved problems of particle physics. On the other hand, neutrino masses 
and other more hypothetical properties such as electromagnetic couplings can play an 
important role in various astrophysical environments. Therefore, astrophysics plays an 
important role at constraining nonstandard neutrino properties. 

Of course, the most exciting recent development is the overwhelming evidence for neu- 
trino oscillations from the atmospheric neutrino anomaly, and notably the zenith-angle 
variation observed in the SuperKamiokande experiment. Besides the near-certainty that 
the phenomenon of neutrino oscillations is real, this high-statistics experiment has also 
opened a new era in neutrino astronomy. It may not be too long until large-scale neutrino 
telescopes observe novel astrophysical sources in the "light" of neutrino radiation. 

In these lectures we cover a number of topics in the area of neutrino astrophysics 
and cosmology which are of current interest to an audience of students who have not 
had much prior exposure to either neutrino physics, astrophysics, or cosmology. At the 
summer school, the lectures were presented on a chalk board, with only a small number of 
viewgraph projections, severely limiting the amount of material that could be covered in a 
few hours. Some of the material was treated in two exercise sessions; some of the exercises 
are integrated into the present notes. Still, these lectures are rather incomplete and give 
only a first impression of the field. 

For a more complete coverage the reader is referred to the excellent textbook by 
Schmitz m, which unfortunately is available only in German. Many of the stellar-evolution 
topics are covered in "Stars as Laboratories for Fundamental Physics" by Raffelt ||. For 
the cosmological questions, the best textbook reference remains the classic by Kolb and 
Turner H. A good overview of many of the relevant issues is provided in two recent text- 
books on astroparticle physics @, [J. Finally, we mention a few recent review articles which 
may be of help to access the field in more depth || [?J §, 3 . 

Chapter 1 treats the production of neutrinos in normal stars, especially in the Sun, 
but leaving out supernova physics — there simply was not enough time to treat this com- 
plicated topic in the lectures. Chapter 2 discusses neutrino oscillations in vacuum and in 
matter. In Chapter 3, experimental strategies are reviewed and some experiments are de- 
scribed in more detail. Chapter 4 is devoted to the connection of neutrinos and cosmology. 
Conclusions are presented in Chpater 5. 
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1 Neutrinos in Normal Stars 



1.1 The Sun 

Ever since it became clear that stars are powered by nuclear fusion reactions and that neu- 
trinos are produced in nuclear reactions, it was also clear that stars are powerful neutrino 
sources. Stellar evolution proceeds through many distinct evolutionary phases [[K], p"T| . 
Stars spend most of their lives on the initial "main-sequence," the Sun is an example, 
where energy is gained from the fusion of hydrogen to helium, i.e. by the net reaction 

Ap + 2e _ -> 4 He + 2u e + 26.73 MeV. (1) 

The detailed reaction chains and cycles depend on the stellar mass which, in turn, influences 
the equilibrium temperature in the interior. In the case of a low-mass star like the Sun, 
hydrogen burning proceeds primarily through the pp chains. The CNO cycle (Bethe- 
Weizsacker cycle), which dominates in more massive stars, contributes only about 1.6%. 

From the particle physics perspective, the solar neutrino flux is perhaps the most im- 
portant example because it has been measured in several different experiments, giving rise 
to the "solar neutrino problem" and thus provides evidence for neutrino oscillations. From 
Eq. ([!]) one sees that two electron neutrinos are produced for every 26.7 MeV of liberated 
nuclear energy. Assuming that the neutrinos themselves carry away only a small fraction 
of this energy, the total solar flux at Earth can be estimated as 

S „ 8.5 x 10 11 MeV cm" 2 s" 1 



6.4 x 10 10 cm" 2 s~\ (2) 



26.7 MeV 26.7 MeV 

where S is the solar constant. 

However, for the interpretation of the experiments, the detailed spectral characteristics 
of the solar neutrino flux are of great importance. In the pp chains, electron neutrinos are 



Table 1: Solar neutrino production in the pp chains. 



Name 


Reaction 


(K) 




Fractional 






[MeV] 


[MeV] 


solar flux 


pp 


p + p^D + e + + v e 


0.26 


0.42 


0.909 


pep 


p + e~+p^D + v e 


1.44 




2 x 10" 3 


hep 


3 He + p -> 4 He + e + + v e 


9.62 


18.77 


2 x 10" 8 


7 Be 


7 Be + e~ -> 7 Li + v e 


(90%) 0.86 




0.074 






(10%) 0.38 






8 B 


8 B -> 8 Be* + e+ + v e 


6.71 


« 15 


8.6 x 10" 5 
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1 Neutrinos in Normal Stars 





Figure 2: Solar neutrino spectrum and thresholds of solar neutrino experiments as indicated 
above the figure (taken from [http:/ /www.sns.ias.edu/^jnb/|) . 



1.2 Photon Dispersion in Stellar Plasmas 3 

produced in six different reactions, giving rise to as many different spectral contributions 
(Table [l], Figs, [l] and |2|), i.e. three monochromatic lines from electron capture reactions and 
three continuous beta-spectra. The sum of the fractional solar neutrino fluxes in Table [l] 
is less than unity due to the small CNO contribution. 

The detailed contribution of each reaction is based on a standard solar model (SSM) fll2f 
which describes the Sun on the basis of a well-defined set of input assumptions. There is 
broad consensus in the literature on the properties of a SSM. However, if the real Sun is 
indeed well represented by a SSM is not a trivial question. The enormous recent progress 
in the field of helioseismology, however, appears to confirm many detailed properties of 
the SSM. One measures the frequencies of the solar pressure modes (p modes) by the 
Doppler shift of spectral lines. One thus probes indirectly the speed of sound in the Sun at 
various depths, i.e. one can reconstruct a sound-speed profile of the Sun which is extremely 
sensitive to the temperature, density and composition profile. 



1.2 Photon Dispersion in Stellar Plasmas 

The nuclear reactions discussed thus far produce electron neutrinos in beta reactions where 
a v e appears together with an e + or at the expense of the absorption of an e~ . However, 
stellar plasmas emit neutrinos also by a variety of processes where vv pairs of all flavors 
appear by effective neutral-current reactions. The most important cases are 

Photo Process 7 + e~ — > e~ + v + u, 

Pair Annihilation e + + — > v + u, 

Bremsstrahlung e~ + (A, Z) — > (A, Z) + e~ + v + u, 

Plasma Process 7 — > v + v. (3) 



The relative importance of these processes depends on the temperature and density of 
the plasma. The energy of the neutrinos produced in these reactions is of the order of 
the temperature of the plasma, in contrast with the nuclear reactions where the neutrino 
energy is determined by nuclear binding energies. In the Sun, the central temperature is 
about 1.3 keV so that thermal neutrino energies are much smaller than those produced in 
the pp chains. The total energy emitted by these processes in the Sun is negligibly small. 

However, in later evolutionary phases, neutrinos produced by plasma processes become 
much more important than nuclear processes. In particular, the plasma process ("pho- 
ton decay") is the dominant neutrino-producing reaction in white dwarfs or the cores of 
horizontal-branch stars or low-mass red giants. This process is noteworthy because it is 
not possible in vacuum due to energy-momentum conservation. In a plasma, on the other 
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hand, the photon acquires a nontrivial dispersion relation ("effective mass") so that its 
decay into neutrinos is kinematically possible. 

We will see that medium-induced modifications of particle dispersion relations are not 
only important for the plasma process, but also for neutrino oscillations in the Sun or 
in other environments. One usually defines a refractive index n re f r which relates wave 
number and frequency of a particle by k = n ve ^u. Refraction in a medium arises from 
the interference of the incoming wave with the scattered waves in the forward direction. 
Therefore, the refractive index is given in terms of the forward-scattering amplitude fo by 

2tt 

n refr = 1 + — nf (u), (4) 

where n is the number density of the scattering targets. This formula applies to any particle 
propagating in a medium, except that fo must be calculated according to the interactions 
of that particle with the medium constituents. 

Photons interact electromagnetically; the dominant scattering process is Compton scat- 
tering on electrons 7 + e~ — > e~ + 7. In the low-energy limit (Thomson scattering) one 
finds /o = —a/m e with a = 1/137 the fine-structure constant. It is then trivial to show 
that the refractive index corresponds to a dispersion relation 

m P 



J 1 - k 2 = ul = — — n e , (5) 



where uj v is the so-called plasma frequency. In the Sun, for example, one finds uj v ~ 0.3 keV 
while in the core of low-mass red giants it is about 9 keV (Exercise |1.4.1| ). The plasma 



frequency plays the role of an effective photon mass. We stress, however, that the general 
photon dispersion relation or that of other particles can not be written in this simple form, 
i.e. in general the effect of dispersion can not be represented by an effective in-medium 
mass. 

1.3 Neutrino Refraction in Media 

We next turn to the neutrino dispersion relation in media. Usually we will be con- 
cerned with very low energies. Therefore, we may take the low-energy limit of the weak- 
interaction Hamiltonian where the propagators for the massive gauge bosons are expanded 
as 

In this limit one obtains the usual current-current Hamiltonian for the neutrino-fermion 
interaction, 

H int = % *Pn»(C v - C AT5 )^7 M (1 - 75)^, (7) 



v r>2 _ i\/r2 - 1\J2 ■ W 
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where G F is the Fermi constant. 

One may then proceed to calculate the dispersion relation in a medium on the basis of 
Eq. (Q). However, in the special case of a current-current interaction the neutrino energy 
shift in a medium can be calculated in a much simpler way. To this end we evaluate the 
expectation value of the current of the background fermions, (V>/7 M (CV — C*a7s The 
axial part (the term proportional to Ca) vanishes if the medium is unpolarized so that 
there is no preferred spin direction. The vector part is equivalent to (jif — nfju^, where n/ 
and nj are the particle and antiparticle densities, respectively, and u is the medium's four- 
velocity. Furthermore, we are only concerned with left-handed neutrino fields for which 
Ihtyu — —^v or (1 — 75)^ = 2ij) v . In summary, the interaction Hamiltonian of Eq. (0) 
amounts to 

V2G F C v (n f - n/K */V7 M Vv (8) 

In the rest frame of the medium we have no preferred direction (no bulk flows) so that 
u = (1,0,0,0). Therefore, a left-handed neutrino in a background medium feels a weak- 
interaction potential 

V = ±V2G F C' v (n f - nj), (9) 
where the lower sign refers to anti-neutrinos. The dispersion relation is 

cu = V + k, (10) 

which evidently does not resemble the one for a massive particle. Therefore, one can not 
define an "effective neutrino mass" in the medium. 

The relevant coupling constants C' v for various background particles are given in Ta- 
ble 0. For most cases, C' v is identical with the neutral-current coupling Cy. However, if / 
is the charged lepton belonging to the neutrino, an additional term with Cy = 1 from the 
Fierz-transformed charged current occurs. For / = v we have a factor 2 for the exchange 
amplitude of two identical particles. 

In an electrically neutral, normal medium we have as many protons as electrons, at 
least if the temperature is low enough that muons and pions are not present, so that 

t/ _i_ /o/-» \-\n n + n e for u e , 

V = ±V2G F x < i , , (11) 

The plus sign is for neutrinos, the minus sign for antineutrinos. For n e < | n n (e.g. in a 
neutron star) the potential produced by the medium is negative for neutrinos and positive 
for antineutrinos. It is important to note that the extra contribution for the electron flavor 
stems from the charged-current interaction in v e + e — > v e + e, which is not possible for 
the other flavors. 
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Table 2: Effective coupling constants for refraction of neutrinos in a medium of background 
fermions. Note that sin 2 Qw — 0.226. 



Fermion / 


Neutrino 


C'y 


Electron 




+| + 2sin 2 9 w 






-\ + 2sin 2 6 H / 


Proton 


V e , Vp, V T 


+\ - 2sin 2 e w 


Neutron 


V e , Up, V T 


i 

2 


Neutrino u a 


v a 


+ 1 




Vb^a 





1.4 Exercises 

1.4.1 Constraints on Neutrino Dipole Moments 

Neutrinos may have anomalous electric and magnetic dipole and transition moments, which 
are small in the standard model but can be large in certain extensions so that they have 
to be constrained. The Lagrangian is 

Ant = ^ Y. (/ i a6^a<VV'& + iivl&k) F ^ i ( 12 ) 

a,b 

where the indices a and b run over the neutrino families, /i a & is a magnetic, e a b an electric 
transition moment, respectively, and a static magnetic or electric dipole moment for a = b. 
(Note that electric dipole moments are CP violating.) These moments are measured in 
units of the Bohr magneton \Ib — e/2m e . F^ u is the electromagnetic field tensor and 

a) Calculate the decay width 7 — >• vv for a neutrino family with a magnetic dipole 
moment //, when photons have an effective mass Wp, as seen in Eq. ([5]). 

b) Calculate the energy loss rate in neutrinos of a nonrelativistic plasma at the temper- 
ature T. 

c) The cores of low-mass red giant stars (about 0.5 M & , solar mass M Q = 2 x 10 33 g) 
have an average density of approximately 2 x 10 5 g cm -3 and are almost isothermal 
at 10 8 K. In order not to delay helium ignition too much, the neutrino loss rate e is 
not allowed to exceed about 10 erg g -1 s _1 . Which limit is obtained for /i? 
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d) This limit is also valid for transition moments, which can lead to decays like v-i — > v\ + 
7. Why does the direct search for these radiative decays of massive neutrinos make 
no sense, provided one believes in the bound for \i obtained above? 

e) Estimate or calculate a similar bound for a hypothetical electric charge (millicharge) 
of a neutrino. 



Hints 

Work in the rest frame of the medium. Show that the squared and spin summed matrix 
element is of the form 

\M | 2 = M a pP a P p , (13) 
where P and P are the four-momenta of the neutrino and antineutrino, respectively, and 

M al3 = V {2K a K p - 2K 2 e* a e p - K 2 g a/3 ) , (14) 

where K = (u>, k) is the photon four momentum and e its polarization four vector. We 
have e*e Q = —1 and e a K a = 0. For the neutrino phase-space integration use Lenard's 
formula 

/ 2^ l| ^ ^-P-P) = ^ (*V* + 2K°K'). (15) 
The decay width is finally found to be 

^ (^ 2 - fc 2 ) 2 

r 7 _> „ = (16) 

Note that the decay would be impossible in vacuum where uo = k. In the present situation, 
however, we may insert the dispersion relation u 2 — k 2 = uj 2 to obtain the decay rate. 

Since every photon decay liberates the energy u in the form of neutrinos, the energy-loss 
rate per unit volume is 

Q = g J pjjAwr^, (17) 

with the photon distribution function (Bose-Einstein) and g = 2 the number of polar- 
ization states. Useful integrals for this exercise are given in Table |A.1| . 

The matrix element and the width for the radiative neutrino decay u 2 — > V\ + 7 (tran- 
sition dipole moment //i 2 ) are 

\Mf = ^UK-P X \K-P 2 ), 

2/2 2\ 3 ( 18 ) 

_ ^ / m 2 -m 1 \ 

where Pj denotes the momentum of neutrino v\ with mass rrii and ii' the photon momentum. 
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1.4.2 Supernova Neutrinos and Refraction 

A type II supernova explosion is actually the implosion of the burnt-out iron core (mass 
around 1.4 M ) of a massive star. This collapse leads to a compact object with nuclear 
density (3 x 10 14 g cm -3 ) and a radius of about 12 km. 

a) What is the gravitational binding energy? 

b) 99% of this energy is emitted in z/s and i>'s of all flavors. When the time for this 
process is 10 s, what is the luminosity (in erg s _1 ) in one neutrino degree of freedom? 

c) The average energy of the emitted neutrinos is 



I 20 MeV other. 

What is the number flux at the time of emission? What is the local neutrino density 
(per flavor) as a function of the radius above the neutron star surface? 

d) At the surface of the neutron star the matter density falls off steeply. Assume that 
is follows a power law p = pn(R/r) p with p = 3 — 7 and pn = 10 14 g cm' 3 . How 
does the electron density compare with the neutrino density during their emission? 
Assume that the medium has as many protons as neutrons. 

e) Compare the weak potential produced by the neutrinos with the one produced by 
normal matter. Assume that the energy flux is the same in all flavors, but the average 
energy is not, see Eq. (|19|). Therefore, only for the electron flavor a difference between 
particles and antiparticles exists and thus a net contribution to the weak potential. 
Another important point is that the v's are emitted almost collinear so that the 
background medium of neutrinos is not isotropic relative to a test neutrino; for an 
exactly pointlike source there would be no contribution at all. 

1.4.3 Neutrino Refraction in the Early Universe 

The "normal" neutrino refractive index is calculated on the basis of the Fermi interaction 
(current-current interaction). It can be interpreted as a weak potential for the neutrinos. 
The medium in the early universe is almost CP-invariant — all particles have almost the 
same number density as their antiparticles. Thus the refractive index nearly cancels to this 
order. A weak potential arises only from the matter-antimatter asymmetry of i] ~ 3 x 10~ 10 
baryons per photon. 



= < 



' 10 MeV 
14 MeV 



for v, 
for v. 



(19) 
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a) Using this value for r], estimate the "normal" refractive index of neutrinos in the 
radiation dominated era before e + e~ annihilation. Use dimensional arguments to 
express n re f r as a function of the cosmic temperature T. (Hint: the number density 
of relativistic thermal particles is proportional to T 3 .) 

b) Which Feynman diagrams contribute to neutrino forward scattering and thus to the 
refractive index? 

c) The gauge boson propagator can be expanded if the momentum transfer Q is small 
relative to the gauge boson mass, 

Ivl z,w Ivl z,w 

The first term provides the Fermi theory of weak interactions. For which diagram 
is the current-current term the exact result? For which diagrams does one have to 
take higher terms into account? 

d) Can you imagine other corrections which might be as important as the propagator 
expansion? 

e) Estimate, again in form of a dimensional analysis, the contribution of the higher 
terms in the early universe. Compare with a). Interpretation? 



Remark 

If the medium consists of neutrinos of all flavors and of e + e~, an exact calculation for 
the CP asymmetric contribution in the early universe gives n rc f r — 1 = ^G 2 F T A /a with 
£ = H 7r(3 — sm2 ®w) srn2 @w — 0.61 for v e or 9 e , while for the other flavors one finds 
f = i|7r(l - sin 2 W ) sin 2 O w ~ 0.17. 



1.4.4 The Sun as a Neutrino Lens 

The neutrino refractive index in media is important for oscillation phenomena. Can it be 
responsible for conventional refractive effects? Estimate the deflection angle of a neutrino 
beam when it hits a spherical body with a given impact parameter. Give a crude numerical 
value for the Sun. Compare with gravitational deflection. 



Hints 

Assume parallel layers of the medium and a beam which propagates at an angle a relative 
to the density gradient. The refraction law informs us that n rc f r sin a = const, where n re f r is 
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the refractive index. Differentiating and some manipulations give for the beam deflection 



where s is a coordinate along the beam and Vj_ the gradient transverse to the local beam 
direction. Since |n re f r — 1| <C 1 one can take n re f r = 1 in the denominator. The deflection 
is so small that to lowest order the beam travels on a straight line. 

The neutrino refractive index is n re f r = l—ml/2E^ — V/E, where V is the weak potential 
of Eq. ([□]). Numerically one finds \p2G F pjm u = 0.762 x 10~ 13 eVp/(g cm -3 ) with the 
mass density p and the atomic mass unit m u . The density at the center of the Sun is about 
150 g cm -3 and the radius of the Sun is R & = 6.96 x 10 10 cm. 

Gravitation affects the beam deflection through the "refractive index" n re f r = 1 — 2$, 
where $ is the gravitational potential. In natural units, Newton's constant is Gn = l/m^ 
with the Planck mass m P i = 1.22 x 10 19 GeV. 



da 
ds 



72 re f r 



(21) 
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2 Neutrino Oscillations 

2.1 Vacuum Oscillations 

If neutrinos have nonzero masses — and thus have properties beyond the standard model — 
they can oscillate between the flavor eigenstates. A flavor eigenstate is operationally defined 
as a neutrino state which appears in association with a given charged lepton. For example, 
the anti-neutrino emerging from neutron decay n —>■ p + e~ + u e is by definition an electron 
anti-neutrino. Likewise, the reaction ^ + p — > n + produces a muon neutrino. Within 
the standard model, flavor-changing neutral currents do not exist, i.e. in a scattering 
process of the form v + n —>■ n + v the out-going neutrino has exactly the same flavor as 
the incoming one. However, in analogy to the quark sector, the flavor eigenstates need not 
be eigenstates of the mass operator. If neutrinos have masses at all, it is generally assumed 
that the mass operator violates the conservation of individual lepton-flavor numbers. 

The simplest example is that of two lepton families. The mass eigenstates Uj, j = 1, 2 
are connected to the flavor eigenstates u a and vp via 

v a \ ( cos 9 sin 9 \ I V\ 
up J I — sin 9 cos 9 J \ u 2 

Depending on the context, Uj may stand for the neutrino field operator or simply for a 
neutrino state, in which case Uj stands for \vj). The mixing matrix is unitary and therefore 
has one nontrivial free parameter, the mixing angle 9. For three (n) families, the mixing 
matrix would have three [\n{n — 1)] nontrivial angles 9j. In addition, for Dirac neutrinos 
it has one [\{n — 2)(n — 1)] CP-violating phase(s), and three [\n(n — 1)] nontrivial phases 
for Majorana neutrinos. However, it can be shown that in oscillation experiments only the 
number of phases given by the Dirac case, i.e. one [\{n — 2)(n — 1)] can be measured. We 
have therefore the same structure as for the Cabbibo-Kobayashi-Maskawa (CKM) matrix 
in the quark sector of the standard model. 

It is easy to see how neutrino oscillations arise if we imagine that a neutrino with 
energy E of a given flavor, for example z/ e , is produced at some location x = 0. It can be 
decomposed into mass eigenstates according to Eq. (|22|) so that 

i^(O) = v e = cos 9 ui + sin 6* u 2 . (23) 

We imagine the mass eigenstates to propagate as plane waves so that each of them evolves 
along the beam as 

v . e -i(Et- Pr x) ( 24 ) 
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with j = 1 or 2. Here 



Pj - = |p,| = (£ 2 - m 2 ) 1 / 2 -E-^, (25) 

where the last approximation holds in the relativistic limit rrij <C E. This is surely justified 
since one expects rrij to be smaller than a few eV and typical energies start at a few MeV. 
Therefore, the original state of Eq. (^) evolves as 

where x — |x|. Next, we can invert Eq. (p2|) to express the mass eigenstates in terms of the 
flavor eigenstates, and insert these expressions into Eq. fl26|). Assuming the other flavor is 
Vfj,, one then finds immediately that the amplitude of the beam evolves as 

(^ e ) = I sin 2 ^ (e^V2E) X _ e i{ml/2E)^ e -iE(t-x) i (2?) 

Squaring this amplitude gives us the probability for observing a ^ at a distance x from 
the production site 

P(^ e -^) = sin 2 2# sin 2 (^^, (28) 

where 5m 2 = m 2 , — m\. A mono-energetic neutrino beam thus oscillates with amplitude 
sin 2 29 and wave number k osc = 5m 2 /4E (Fig. |^). The maximum effect occurs for 9 = it /A. 
One usually defines the oscillation length 

4nE E 1 eV 2 

Lsc = n = 2.48 m -. 29 

Sm 2 1 MeV Sm 2 V ; 

The neutrino beam has returned to its original state after traveling a distance £ osc . The 
probability for finding the neutrino in its original state after traveling a distance x is 
P(i/ e ->i/ e ) = l-P(i/ e ->!/„). 

Note that oscillations would be impossible for completely degenerate masses (m 2 = m 2 ), 
including the case of vanishing neutrino masses, or for a vanishing mixing angle. 

The generalization of these results to three or more families is straightforward but 
complicated; it can be found in many textbooks, e.g. in Ref. [p]]. Equation (|28"D then reads 



j>i 



1 - exp ( -i^ x 



(30) 



with Aij = mf — m 2 . In general P(a —>■ (3) ^ P(f3 —>■ a) due to the complex phase of 
the unitary mixing matrix U, offering a possibility to measure this phase. A convenient 



2.2 Oscillation in Matter 
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prob(i/ e ^i/^) 




Figure 3: Oscillation pattern for two-flavor oscillations (neutrino energy u>, distance z). 



parametrisation of U is 



/ 



U 



Cl 



sic 3 



-S1C2 C1C2C3 - s 2 S3e 



C1C2S3 + s 2 c 3 e l 



(31) 



V sis 2 -C1S2C3 - C2S 3 e l<5 -C1S2S3 + c 2 c 3 e lS J 



with Si - 
diag(e* Al . 



sin^j and q = cos^j. For Majorana neutrinos U has to be multiplied with 



Thus far we have assumed that the neutrinos are mono-energetic and the sources and 
detectors are pointlike. Since nature is not so kind as to provide us with these simple cases, 
one has to convolute these formulas with energy and distance distributions. Naturally, these 
effects smear out the signature of oscillations, thereby complicating the interpretation of 
the experiments. A given experiment or observation usually provides exclusion or evidence 
regions in the parameter plane of sin 2 29 and 5m 2 . 



2.2 Oscillation in Matter 
2.2.1 Homogeneous Medium 

Neutrino oscillations arise over macroscopic distances because the momentum difference 
Sm 2 /2E between two neutrino mass eigenstates of energy E is very small for neutrino 
masses in the eV range or below and for energies in the MeV range or above. Wolfenstein 
was the first to recognize that the neutrino refractive effect caused by the presence of 
a medium can cause a momentum difference of the same general magnitude, implying 
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that the extremely small weak potential for neutrinos in a medium can modify neutrino 
oscillations in observable ways. 

In order to understand how the neutrino potential enters the oscillation problem, it is 
useful to back up and derive a more formal equation for the evolution of a neutrino beam. 
To this end we begin with the Klein-Gordon equation for the neutrino fields 

(d 2 - V 2 + M 2 )^ = (32) 

where in the general three-flavor case 



M 2 



fm\ 

ml j and * = | v 2 | . (33) 
V m 2 . 




If we imagine neutrinos to be produced with a fixed energy E at some source, their wave 
functions vary as e~ %Et so that their spatial propagation is governed by the equation 

(-£ 2 -<9 2 + M 2 )^ = 0, (34) 

where we have reduced the spatial variation to one dimension, i.e. we consider plane waves. 

In the relativistic limit E ^> m 2 we may linearize this wave equation by virtue of 
the decomposition (— E 2 — <9 2 ) = — (E + id x ){E — id x ). Since —id x Vj = pjVj with pj = 
(E 2 — m 2 ) 1 / 2 ~ E it is enough to keep the differential in the difference term, while replacing 
it with E in the sum, leading to (— E 2 — d 2 ) — > —2E(E + id x ). Therefore, in the relativistic 
limit the evolution along the beam is governed by a Schrodinger-type equation 

id x * = (-E + n)*, «=^- (35) 

Usually this sort of equation is written down for the time-variation instead of the spatial 
one so that it looks more like a conventional Schrodinger equation. However, in the problem 
at hand we ask about the variation of the flavor content along a stationary beam, so that 
it is confusing to use a differential equation for the time variation which then has to be 
re-interpreted as describing the evolution along the beam. Either way, the main feature of 
this equation is that it is a complex linear equation involving a "Hamiltonian" matrix Q; 
the term — E contributes a global phase which is irrelevant for the oscillation probability. 

The potential caused by the medium is then easily included by adding it to the Hamil- 
tonian 

M 2 

n - n M = — + v (36) 

where V is a matrix of potentials which is diagonal in the weak-interaction basis with the 
entries given by Eq. flTT|). 
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As a specific example we now consider two-flavor mixing between v e and z/ M , and write 
down the Hamiltonian in the weak interaction basis. It is connected to the mass basis by 
virtue of the unitary transformation v a = U a iVi of Eq. fl2"2"|) . The squared mass matrix then 
transforms as UM 2 U\ leading in the weak basis explicitly to 



9 £ 

M = — 
2 



5m 2 



-cos 26 
sin 26 



sin 26 ' 
cos 26 



(37) 



where £ = m\ + m\ and Sm 2 
"Hamiltonian" 



2 

mf. Including the weak potential then leads to the 



n 



5m 2 ( — cos 26 



M 



+ V2d 



+ 








- 2 77- ri 



+ - sin 26 
2 




(3* 



which governs the oscillations in a medium. The first term which is proportional to the unit 
matrix produces an irrelevant overall phase so that the medium effect on the oscillations 
depends only on the electron density n e , i.e. on the difference between the weak potentials 
for v e and v^. Recall that this difference arises from the charged-current piece in the v e 
interaction with the electrons of the medium. 

The meaning of this complicated-looking expression becomes more transparent if one 
determines the "propagation eigenstates," i.e. the basis where £7m is diagonal. In vacuum 
we have M 2 = 2EQ so that in the medium one may define an effective mass matrix by 
virtue of = 2EQm- The eigenvalues of this matrix are found in the usual way by 



solving det(2EQ 



M 



m 



M 



0. The two roots and their difference are found to be 
= \ (S + 2 v / 2G F (n e - n n )E =p 5m 2 M 



(39) 



5m 2 



1/2 



mi — mi. We 



X1 - [(Sm 2 ) 2 + AEG F n e (2EG F n e - V2 5m 2 cos 26 

In vacuum where n e = n n = these expressions reduce to m\ 2 and 5m 2 
stress that the in-medium effective squared "masses" should not be literally interpreted as 
effective masses as they depend on energy; may even become negative. The "effective 
masses" are simply a way to express the dispersion relation in the medium. 

The transformation between the in-medium propagation eigenstates and the weak- 
interaction eigenstates is effected by a unitary transformation of the form Eq. (|22|) with 
the in-medium mixing angle 

sin 20 

tan2# M = — n r (40) 



cos 26- A 



16 



2 Neutrino Oscillations 



which is equivalent to 

9 „ sin 2 29 

sin 2 29 M = 5—. 41 

(cos 2^- A) 2 + sin 2 20 v ; 

Here 

itB 2 3 ^g g , = 1Jaxl0 - T _^L E eVl 
om z gem - - 3 JVleV om 2 

where F e is the electron number per baryon and p the mass density. 

With these results it is trivial to transcribe the oscillation probability from the previous 
section to this case, 

f xbm\ x 



Evidently we have 



Pu(ve -> ^) = sin 2 29 M sin 2 I I . (43) 



_ 4vr£ _ sm2gM ^ 

*M — T 2~ ~~ • a *-vac, (44 J 

07?% sm 26 1 



for the in-medium oscillation length. 
2.2.2 MSW Effect 



The mixing angle in matter, sin 29m-, is a function of n e and E. It becomes maximal 
(9m = 7t/4) when A = A R = cos 29. Here ^m 2 ^ has a minimum, the oscillation length £ M 
a maximum. Even if the mixing angle in vacuum is small, on resonance the in-medium 
mixing is maximal, independently of the mixing angle in vacuum. 

The most interesting situation arises when a neutrino beam passes through a medium 
with variable density, the main example being solar neutrinos which are produced at high 
density in the solar core. Considering the case m\ ~ and 9 small, we find for high density 
(A > A R ) that 9 M — \, implying 

v X m - -Vfi, m\ M ~ 0, 

z/ 2M ^ +^ e , rn\ M - A 5m 2 . (45) 

On the other hand, for low densities (A <C A R ) we have vacuum mixing {9m — 9 <C 1) so 
that 

v lM ~ +v e , m\ M ~ 0, 

z/ 2M ^ rnl M ^ rn\. (46) 

Therefore, the propagation eigenstate v 2 m which at high density is approximately a u e 
turns into a at low density, and vice versa. Therefore, if the neutrino is born as a v e , 
and if the density variation is adiabatic so that the neutrino can be thought of being in a 
propagation eigenstate all along its trajectory, it will emerge as a v^. 



2.3 Exercises 
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In order for this resonant conversion to occur, two conditions must be met. First, the 
production and detection must occur on opposite sides of a layer with the resonant density. 
In the Sun, the neutrinos are produced at high density so that we need to require 

A(place of production) > A R = cos 29. (47) 

This can be rewritten as a constraint on the neutrino energy, 

_ „ 5m 2 cos 29 K , , , r „ 5m 2 g cm -3 . . 

E>E = — = = 6.6 x 10 6 MeV cos2fl — — ^— . 48 

2V2G F n e eV 2 Y e p v ; 

Second, for the neutrino to stay in the state z/ 2m , the density gradient has to be moderate, 
i.e. the density variation must be small for several oscillation lengths. This condition can 
be expressed as £^ ^ Vlnn e . One often defines the adiabaticity parameter 

_ 5m 2 sin 2 29 1 
7 ~ 2E cos 29 |Vlnn e | ^ ' 

so that adiabaticity is achieved for 7 < 1. This condition must be met along the entire 
trajectory. As the oscillation length is longest on resonance when the mixing is maximum, 
the adiabaticity condition is most restrictive on resonance. In general one finds a triangle 
in sin 2 29-5m 2 space where these conditions are fulfilled — see Exercise [2.3.2| . In Fig. |] the 
triangular contours show the range of masses and mixing angles where the solar v e flux is 
reduced to the measured levels for experiments with different spectral response, i.e. which 
measure different average neutrino energies. 



2.3 Exercises 

2.3.1 Polarization vector and neutrino oscillations 

Neutrino oscillations are frequently described by a Schrodinger equation of the form 

M 2 

= with Q = p H , (50) 

2p 

with p the neutrino momentum, M the mass matrix, and ^ a column vector with two 
or more flavors. For two generations, the relation between flavor and mass eigenstates is 
given by Eq. (p2|). Instead of the state vectors, however, one can work with the 2x2 density 
matrix in flavor space which is defined by 

Pah = %^a, (51) 

where the indices a and b run, for example, over v e and or over 1 and 2 in the mass basis. 
With the help of the density matrix we can find an intuitive geometric interpretation of 
oscillation phenomena. In addition, one can treat statistical mixtures of states, i.e. when 
the neutrinos are not characterized by pure states. 
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a) Show that the equation of motion is: ip = [ft, p] = [M 2 , p]/2p. 

b) Write the mass matrix in the form M 2 /2p = Vq — |V ■ o and show, that in the flavor 
basis 

/ sin 29 \ 



ml + m\ 2tt 
V = —. and V = 

4p ^OSC 





^ cos 29 J 



Atx p 

with u osc = —= j ■ ( 52 ) 

mi — mi 



The vector V is thus rotated against the 3-axis with the angle 29. Has this orientation 
in the 1-2 plain a physical meaning? 

c) Express the density matrix in terms of a polarization vector in form of p = |(l+P-er). 
Physical interpretation of its components? 

d) Which property of P characterizes the "purity" of the state, i.e. when does the density 
matrix describe pure states, when maximally incoherent mixing? 

e) Show that the equation of motion is a precession formula, P = V ■ P. Obtain the 
oscillation probability for an initial u e . 

f) The energy of (non-mixed) relativistic neutrinos in a normal medium is E = p + 
(m 2 /2p) + Vmed- Here V me d is given by Eq. (|ll|). What is P in the medium? What 
is the mixing angle in the medium? 

g) In a medium consisting of neutrinos (supernova, early universe) one can not distin- 
guish between a test neutrino and a background neutrino, so that oscillations with 
medium effects are in general nonlinear. What is the advantage of the density matrix 
formalism in this situation? 



2.3.2 MSW Effect in the Sun 

The conditions for the MSW effect are given by the Eqs. ( |18D and (f£5|). Determine the 
region in sin 2 29-5m 2 space, where one expects almost complete flavor inversion, i.e. the 
MSW triangle. For this purpose, assume that all solar neutrinos are produced with an 
energy of E = 1 MeV and that the density profile of the Sun is approximated by an 
exponential of the form n e = n c exp(— r/Ro), with a scale height of i?o = -R©/10.54 and a 
density at the center of n c = 1.6 x 10 26 cm -3 . 




Figure 4: Limits and evidence for neutrino oscillations (Figure courtesy of Max Tegmark). 
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3 Experimental Oscillation Searches 
3.1 Typical Scales 

We now turn to a discussion of some experimental strategies for the detection of neutrino 
oscillations. The most widely used neutrino sources are the Sun, the Earth's atmosphere 
where neutrinos emerge from cosmic-ray interactions, or man-made devices such as re- 
actors and accelerators. One distinguishes between appearance and disappearance exper- 
iments. In the former, one searches for the appearance of another flavor than has been 
produced in the source, while the latter are only sensitive to a deficit of the original flux. 

From Eq. fl2~8|) one finds that an experiment with typical neutrino energies E and a 
distance L between source and detector is sensitive to a minimal value of the mass-squared 
difference of 

E 

(Sm 2 ) min ~ — . (53) 

Therefore, different experiments probe different regions of the mass sector. In Table [3] we 
give some examples. 



Table 3: Characteristics of typical oscillation experiments. 



Source 


Flavor 


E 


[GeV] 


L [km] 


(<5m 2 )min [eV 2 ] 


Atmosphere 


(-) (-) 


10" 


-l 


10 2 


10. . .10 4 


10~ 6 


Sun 




io- 


-3 


io- 2 


10 8 


io- 11 


Reactor 


Ve 


io- 


-4 


io- 2 


io- 1 


10~ 3 


Accelerator 


(-) 


10" 


-1 


1 


10 2 


io- 1 



3.2 Atmospheric neutrino experiments 

When cosmic rays, i.e. protons and heavier nuclei interact with the Earth's atmosphere 
they produce kaons and pions, which in turn decay into muons, electrons and neutrinos. 
Since the initial state is positively charged one has more neutrinos than anti-neutrinos, 
but the experiments are insensitive to this effect. On the other hand, the flavor can be 
very well measured; from the simple production mechanism one expects 



N(^) : N(u e ) = 2:1. 



(54) 



3.2 Atmospheric neutrino experiments 
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This ratio depends on the energy of the measured neutrinos since the lifetime of high- 
energy muons is increased by their Lorentz factor so that they may hit the ground before 
decaying. One often uses the double ratio 



R 



(N(n)/N(e)) t 
(NM/N{e)) 



(55) 



MC 



meaning the ratio of the measured flavor ratio divided by the expectation from Monte- 
Carlo simulations. The double ratio has the advantage that the uncertainty of the overall 
absolute flux cancels out. The average neutrino energy is found to be (E„) ~ 0.6 GeV. 

A textbook example of an experiment of this kind is SuperKamiokande (SK) an 
underground detector consisting of 50 kt of water, surrounded by 11000 photomultipliers. 
Neutrinos react with the protons and neutrons of the target and produce electrons and 
muons. These charged particles are identified by their cones of Cherenkov light which are 
fuzzier for the electrons. Since cosmic rays are distributed almost isotropically and the 
atmosphere is spherically symmetric, one expects the flux to be the same for down or up- 
coming neutrinos. However, it was found that up-coming muon neutrinos are significantly 
suppressed. This up-down asymmetry is shown in Fig. |5|. Plotted is the momentum of the 




1 10 

Momentum (GeV/c) 



Figure 5: Up-down asymmetry at SK. The hatched region is the expectation without oscil- 
lation, the dots the measurements, while the dashed line represents the best-fit oscillation 
case. (Figure from JT3].) 
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charged leptons against 

where U and Z) are the number of upward and downward going events, i.e. their zenith 
angle is larger or smaller than about 78°, respectively. The dashed line corresponds to the 
best-fit mixing parameters of SK, 5m 2 = 2.2 x 10~ 3 eV 2 and sin 2 29 = 1. Roughly the 
same parameters are found by similar experiments, like 1MB |TJJ, MACRO [ITH] or Soudan 



Tfi| . For the double ratio R values between 0.4 and 0.7 were found. 

Therefore, the z^'s probably oscillate into z/j-'s. The other possibility, oscillation into 
sterile neutrinos v s is disfavored because the observed rate of the NC process vN — ► vk°X 
is about as expected. Sterile neutrinos by definition do not take part in such reactions. 
In addition, other properties such as the energy distributions of the final-state charged 
leptons tend to confirm the v^—Vt interpretation. 

3.3 Accelerator Experiments 

Evidence for oscillations were present before the SuperKamiokande results. The LSND 



17fl collaboration used a 800 MeV proton beam colliding on a water target so that pions 
were produced. The n~ were captured in a copper block while the n + decayed into /i + z/^. 
These in turn decayed into a positron and two neutrinos. Therefore, one expects the same 
number of u^, and v e , but no u e . In a scintillation detector 30 m behind the source 
they looked for u e in the reaction v e + p — > e + + n. The experimental signature is thus the 
Cherenkov cone from the positron and a photon from the reaction 'j+p — > D + 7 (2.2 MeV). 

The LSND collaboration measured an excess of about 40 events above the background. 
The interpretation as oscillations, however, is controversial because the very similar KAR- 
MEN [^] experiment sees no events in about the same parameter range. On the other 
hand, KARMEN will not be able to exclude the LSND results. The remaining parameter 
range where a consistent interpretation as oscillations remains possible is shown in Fig. [|, 
i.e. 5m 2 ~ 1 eV 2 and sin 2 26 ~ 10~ 2 . 



3.4 Reactor Experiments 

In power reactors, nuclear fission produces v e with energies of typically several MeV. These 
energies are too low to produce a charged mu or tau lepton in the detector so that reactor 
experiments are always disappearance experiments. The detection reaction is v e + p — > 
e + + n. 

Thus far, none of the reactor experiment gives evidence for oscillations. However, they 
have produced very important exclusion areas in the sin 2 29-5m 2 space. Most importantly, 




Figure 6: Exclusion plot of the CHOOZ experiment; the area to the right of the lines is 
excluded. Also shown is the allowed parameter space of Kamiokande. The SK has shifted 
these values towards lower masses, yet they still lie inside the forbidden area. (Figure from 



the CHOOZ experiment [|19[] has excluded a large range of mixing parameters (Fig. |6[) so 
that a putative v^-v^ oscillation interpration of the SuperKamiokande results is inconsistent 
with their limits. 



3.5 Solar Neutrino Experiments 

As discussed in Chapter 1.1, there are six different solar neutrino reactions in the pp 
chain with six different energy spectra. Different experiments measure neutrinos from 
different reactions because they have different energy thresholds and different spectral 



response characteristics. Some of the experiments are Homestake pC§ which uses the 
detection reaction v e + 37 C1 — > 37 At + e~ (threshold 814 keV), the gallium experiments 
GALLEX P], and SAGE Pj which use v e + 71 Ga -»• 71 Ge + e~ (threshold 232 keV) and 
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( Super) Kamiokande |23], Q where the elastic scattering on electrons v e + e~ — > u e + e~ 
is used (threshold 5 MeV). These experiments have in common that they are deep under 
the Earth (typically some 1000 m water equivalent) to eliminate backgrounds from cosmic 
radiation. 

Soon after the first experiments were started it was found that one half to two thirds 
of the neutrino flux was missing. This "solar neutrino problem" is illustrated by Fig. |7] 
where the measured rates of the chlorine, gallium, and water Cherenkov experiments are 
juxtaposed with the predictions in the absence of oscillations. It is important to note that 
the flux suppression is not the same factor in all experiments. Rather, it looks as if there 
was a distinct spectral dependence of the neutrino deficit. 



Total Rates: Standard Model vs. Experiment 
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Figure 7: Solar neutrino measurements vs. theoretical flux predictions. (Figure taken from 
http:/ /www.sns.ias.edu/~jnb/| .) 
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In particular, it seems as if the 7 Be neutrinos do not reach Earth. Best-fit solutions for 
SSM flux variations even yield negative 7 Be fluxes unless neutrino oscillations are taken into 
account. Possible non-oscillation explanations seem to be unable to explain this scenario 
in the light of helioseismology which shows excellent agreement with the SSM. The fact 
that the nuclear cross sections necessary for the SSM are not known for the relevant (low) 
energies but have to be extrapolated from higher energy experiments also can not give the 
measured rates. Temperatures different from the SSM assumption are constrained by the 
crucial T dependence of the fluxes, namely <p u {B 8 ) oc T 18 and ^(Be 7 ) oc T 8 which do not 
allow to reduce the beryllium flux by a larger factor than the boron flux. 

The solar neutrino measurements can be consistently interpreted in terms of two-flavor 
oscillations. In contrast to the atmospheric and reactor results, several different solutions 



exist as shown in Fig. f|. Typical best-fit points are [25|] 



{Sm 2 (eV 2 ), sin 2 29) 



' (5.4 x 1(T 6 , 6.0 x 10- 3 ) SAMSW 

(1.8 x 10~ 5 , 0.76) LAMSW ( 57 ) 
(8.0 x 10" 11 , 0.75) VO 



Here, SAMSW and LAMSW denote the small-angle and large-angle MSW solution, re- 
spectively, while VO refers to vacuum oscillation. 

Strategies to decide between these solutions include more precise investigation of the 
electron energy spectrum (SAMSW/LAMSW) or seasonal variations (VO). A comparison 
of NC and CC events in the SNO experiment |26| will clarify if the v e oscillate into sequential 
or sterile neutrinos. 



3.6 Summary of Experimental Results 

With three different masses there can be only two independent Sm 2 values. Should all 
experiments with evidence for oscillations be confirmed then there must be a fourth neu- 
trino, which has to be sterile because LEP measured the number of sequential neutrinos 
with m u < 45 GeV to be 3. The current results can be summarized as 

(10- 3 , > 0.7) 

(10~ 5 , 10- 3 ) 
(10- 5 ,> 0.7) 
(10- 10 ,>0.7) 
(1,10-3) 
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Many authors tend to ignore the LSND results to avoid the seemingly unnatural possibility 
of a sterile neutrino. On the other hand, if it were found to exist, after all, this would be 
the most important discovery in neutrino-oscillation physics. 

How to choose the masses to get the observed differences is a topic of its own. The most 
interesting aspect of the mass and mixing scheme is the influence on neutrinoless double 
beta decay. We refer to [[| for a review of that interesting issue. For later need in the next 
chapter it suffices to say that if the atmospheric mass scale is interpreted in a hierarchical 
mass scenario (5m 2 ~ m 2 ) we have one mass eigenstate of about 0.03 eV. 

We close this chapter with a few words on future important experiments. For solar 



physics, besides the SNO experiment, the BOREXINO |27j experiment will measure the 
crucial 7 Be flux to a higher precision than its predecessors. MiniBoone [^8] is designed 
to close the LSND/KARMEN debate and thus confirm or refute the need for a sterile 
neutrino. 

Also planned are so-called long baseline (LBL) accelerator experiments from KEK to 
Kamioka (K2K, H), from Fermilab to Soudan (MINOS, ||) and from CERN to Gran 
Sasso (ICARUS, |31[] ) probing mass ranges down to 5m 2 ~ 10~ 3 eV and testing possible 
CP violation in the neutrino sector. In the years to come, exciting discoveries are to be 
expected. 
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4 Neutrinos in Cosmology 

4.1 Friedmann Equation and Cosmological Basics 

The neutrino plus anti-neutrino density per family in the universe of 113 cm -3 is compara- 
ble to the photon density of 411 cm -3 . Therefore, it is no surprise that neutrinos, especially 
if they have a non-vanishing mass, may be important in cosmology. To appreciate their 
cosmological role, we first need to discuss some basic properties of the universe. 

On average, the universe is homogeneous and isotropic. Its expansion is governed by 
the Friedmann equation, 

H 2 = 8 -^G NP - \, (59) 

where is Newton's constant, p the energy density of the universe, and H = a/a the 
expansion parameter with a(t) the cosmic scale factor with the dimension of a length. The 
present-day value of the expansion parameter is usually called the Hubble constant. It has 
the value 

H = h 100 km s" 1 Mpc" 1 , (60) 

where h = 0.5-0.8 is a dimensionless "fudge factor." 

The constant k determines the spatial geometry of the universe which is Euclidean 
(flat) for k = 0, and positively or negatively curved with radius a for k = ±1, respectively. 
The universe is spatially closed and will recollapse in the future for k — +1. For k = 
or —1 it expands forever and is spatially infinite (open), assuming the simplest topological 
structure. However, even a flat or negatively curved universe can be closed. An example 
for a flat closed geometry is a periodic space, i.e. one with the topology of a torus. 

Equation (|59| ) must be derived from Einstein's field equations. However, it can also be 
heuristically derived by a simple Newtonian argument. Since the universe is assumed to be 
isotropic about every point, we may pick one arbitrary center as the origin of a coordinate 
system. Next, we consider a test mass m at a distance R(t) = a(t) r from the center, and 
assume that the homogeneous gravitating mass density is p. The total energy of the test 
mass is conserved and may be written as E = —^Km. On the other hand, 

^ m 1-9 GmMtti , , 

E = T + V = -mR 2 - — , (61) 
2 ti 

where M = R 3 p47r/3 is the total mass enclosed by a sphere of radius R. With K = kr 2 
the Friedmann equation follows. On the basis of a Galileo transformation one can easily 
show that the Friedmann equation stays the same when transformed to another point, i.e. 
the expansion looks isotropic from every chosen center. 

Some basic characteristics of the expanding universe are easily understood if we study 
the scaling behavior of the energy density p. Nonrelativistic matter ("dust") is simply 
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diluted by the expansion, while the total number of "particles" in a co-moving volume, of 
course, is conserved. Therefore, when matter dominates, we find p oc a~ 3 . 

For radiation (massless particles), the total number in a comoving volume is also con- 
served. In addition, we must take the redshift by the cosmic expansion into account. The 
simplest heuristic derivation is to observe that the cosmic expansion stretches space like a 
rubber-sheet and thus stretches the periodic pattern defined by a wave phenomenon. Thus 
the wavelength A of a particle grows with the cosmic scale factor a, implying that its wave 
number k and thus its momentum scale as a -1 . For radiation we have uj = k so that the 
energy of every quantum of radiation decreases inversely with the cosmic scale factor. In 
summary, the energy density of radiation scales as a~ A . 

For thermal radiation (blackbody radiation) we may employ the Stefan-Boltzmann 
law p oc T 4 to see that T oc a -1 . The temperature of the cosmic microwave background 
radiation is a direct proxy for the inverse of the cosmic scale factor. 

The density in Friedmann's equation, therefore, decreases at least with a -3 so that at 
late times the curvature term takes over if k — ±1. Of course, it is frequently assumed 
that the universe is flat, and certainly the curvature term, if present at all, may only begin 
to be important today. For k — — 1, H 2 will never change sign, so the universe expands 
forever. For k = +1, the expansion turns around when H 2 = 0, i.e. when a~ 2 = (8tt/3)GnP- 
Either way, at very early times (large temperature, small scale factor), radiation dominates. 
Therefore, in the early universe we may neglect the curvature term. 

A "critical" or Euclidean universe is characterized by k = 0. In this case we have 
a unique relationship between p and H . The density corresponding to the present-day 
expansion parameter Hq is called the critical density 

p c = ML = h 2 1.88 x lO" 29 g cm" 3 . (62) 
87rGjv 

The last number is the current value of the critical density, denoted p . It translates to 
about 1CT 5 protons per cm 3 . The density is often expressed as a fraction of the critical 
density by virtue of Q = p/p c . A value of ft — 1 corresponds to a flat universe. 

The cosmic background radiation, which was first observed in the 1960s, is equivalent 
to the radiation of a black body with temperature T = 2.726 K. Therefore, today's number 
density of photons is 



n ~< = 2 / = t3 - 411 - 5 cm " 3 > (63) 

where f(uo) = (e u l T — l) -1 is the Bose-Einstein distribution function. The energy density 
of photons (or in general bosons) in the universe is calculated by 
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where we have introduced the effective number of boson degrees of freedom g-Q for later 
use. Therefore, 

fi 7 = ^ ~ 4.658 x 10~ 34 g cm~ 3 /Po ^ 2.480 x 10~ 5 br 2 (65) 
Po 

is the contribution of microwave photons to the critical density. 

The fraction of baryons relative to the number of photons in the present-day universe 

is 

77 = — = 7710 x 1(T 10 ~ 3 x 10- 10 , (66) 



n 7 



so that the contribution of baryons to the critical density is merely 

n B = = = ^ 3 6 x 10 - 3 h -2 ^ 0M h -2 (67) 

Po Po Po 

where we took the proton mass as an average baryonic energy. 

4.2 Radiation Epoch 

In the hot early universe, neutrinos should have been in thermal equilibrium so that one 
expects a cosmological neutrino sea in analogy to the cosmic microwave background. The 
number density of one massless thermal neutrino generation is 

= 2 /|^' = § T3 < (68) 

where we have to use the Fermi-Dirac distribution f(u) = (e w / T + 1) _1 for fermions with a 
vanishing chemical potential. Equation (6~8") counts two internal degrees of freedom, i.e. it 



counts left-handed (interacting) neutrinos and (right-handed) anti-neutrinos for a given 
family. The energy density is 

with the effective number of fermion degrees of freedom g-p. The effective number of thermal 
degrees of freedom is defined by 

7 

9* = E 9j + g E 9j, (70) 

j'=bosons j'=fermions 

so that the total energy density in radiation is 

Prad = #*— T 4 . (71) 
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Another important quantity is the entropy density in radiation S = (p + p)/T. Since for a 
relativistic gas p = p/3 we have 

2tt 2 

S = g, — T\ (72) 

The energy and entropy densities are simply functions of the temperature and of the particle 
degrees of freedom that are thermally excited at that temperature. 

To get a feeling for g*, we consider some characteristic temperatures. At T < 100 MeV, 
below the QCD phase transition, there are photons, e ± , and three flavors of (anti) neutrinos, 
giving g m = 2 + 7/8 x 10 = 43/4 = 10.75. At T = 300 MeV, significantly above the 
QCD phase transition, in addition muons (4 fermionic degrees of freedom), eight gluons 
(2x8 bosonic degrees of freedom) and up, down, and strange quarks (3x4x3 fermionic 
degrees of freedom, counting 3 color degrees of freedom each) were present, resulting in 
g* = 18 + 7/8 x 44 = 113/2 = 56.5. The number of degrees of freedom drops dramatically 
around the QCD phase transition when quarks and gluons condense to confined states of 
mesons and baryons. 

The Friedmann equation for the radiation epoch in the early universe can be written 
in the form 

47T 3 

H 2 = -rrG N g*T 4 . (73) 
45 

We can now calculate a relationship between cosmic age and temperature in the radiation 
epoch. We use p = p (R /R) 4 = bR~ 4 and H 2 = %n /^G^bR' 4 which leads to 



^) = (f^&) 1/2 i (74) 



and thus to 



R(t)=(^G N b) 1/ \ 1 ' 2 . (75) 



3 



With Eq. (|TTD for p(T) we find 



T = V^rV^ 1.5 6 MeV (— • (76) 



\16ir 3 G N g*J '" V t 

Therefore, the universe had a temperature of about 1 MeV when it was about 1 s old. 



4.3 Present— Day Neutrino Density 

When neutrinos and photons are in thermal equilibrium in the early universe, they are 
both characterized by the same temperature. However, the neutrino temperature today is 
thought to be lower than that of the cosmic microwave background. The reason for the 
difference is that when the temperature has fallen below the electron mass m e (t ~ 10 s) 
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the annihilation of electrons and positrons heats the photon gas by virtue of e + e~ — > 27. 
As we will see later, neutrinos have already decoupled at that time so that the rate for 
e + e~ — ► vv is too slow to be of importance. 

The electron-positron-photon plasma is so tightly coupled that one may assume that it 
is always close to thermal equilibrium throughout the annihilation process, i.e. the disap- 
pearance of the pairs is an adiabatic process. Therefore, the entropy of the gas will 
go over to the photons. If we denote the photon temperature before this process as T\ and 
afterward as T 2 , and if the value of g* for the coupled species before and after annihilation 
is #1,2, respectively, entropy conservation implies 

9l Tf = g 2 Tl (77) 

Before annihilation, the participating species are and photons so that g* = 2 + 4x7/8 = 
11/2, while afterward we have only photons with g* = 2. Before annihilation, we have 
T v = T 7 so that the ratio TijT\ can be interpreted as T 7 /T y after annihilation. Therefore, 
we find the famous result 

T v = (yj) r 7 . (78) 

Thus, the present-day temperature of the neutrino background is predicted to be T v ~ 
1.946 K ~ 1.678 x 10" 4 eV. Naturally, its direct detection is an impossible task. 

The number density of one neutrino familiy is given by Eq. fl6"8"|) . Inserting the above 
numbers yields 

n uP ~ 337.5 cm" 3 (79) 
for all three generations. If all neutrinos are relativistic today, all three families contribute 

P^ = ^(tt) — P 7 ^ 3.174 x IO- 34 g cm" 3 (80) 
8\11/ g 7 

to the cosmic energy density. 

We have assumed that the neutrinos are relativistic, which is the case when their mass 
is smaller than their average energy 

(E v ) = ^~3T,~5x 1(T 4 eV. (81) 

If the indications for neutrino oscillations discussed in the previous section are correct, at 
least one of the mass eigenstates exceeds around 0.03 eV and thus is not relativistic today. 

But even if m v is bigger than 5 x 10 -4 eV, the number density is still given by Eq. ([79]) 
since at the point of e annihilation the neutrinos were relativistic. The energy density 
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today, though, changes for one nonrelativistic flavor to p u p = n uP m u . The contribution to 
the critical density today is then Q uP = p uP /p c or 

^E^lv- (82) 

flavors 

We may use this result to derive the famous cosmological limit on the neutrino masses. It 
turns out that the lower limit on the age of the universe indicated by the age of globular 
clusters implies something like Vth 2 ~ 0.4. Since Vt vD < we have 

m„~37 eV. (83) 

flavors 

If the experimental indications for neutrino oscillations are correct, the mass differences 
are very small. Therefore, a neutrino mass in the neighborhood of this limit would require 
that all flavors have nearly degenerate masses, implying that we may divide this limit by 
the number of flavors to obtain a limit on the individual masses. Therefore, 

m v ~ 12 eV (84) 

applies to all sequential neutrinos. 

If the current indications for oscillations from the atmospheric neutrino anomaly are 
correct so that one mass eigenstate exceeds about 0.03 eV we find that VL u ph ~ 3 x 10" 4 
or, with h ~ 0.8, we have 10 4 , not much less than the luminous matter of the 

universe. 

One may wonder if any of these results change if neutrinos are Dirac particles so that 
there are actually 4 degrees of freedom per flavor. However, the sterile (right-handed) com- 
ponents will not be thermally excited in the relevant epochs of the early universe because 
they are too weakly coupled if the masses are as small as indicated by the cosmological 
limit. It is important to realize that g* is the effective number of thermally excited de- 
grees of freedom, which does not necessarily include all existing particles. For example, 
the massless gravitons are another species which are not thermally excited anywhere near 
the QCD or e ± annihilation epochs so that they, too, have not appeared in our particle 
counting for g*. 



4.4 Big Bang Nucleosynthesis 

How do we know that the cosmological neutrino background actually exists? After all, its 
direct detection is an impossible task with present-day experimental means. However, the 
cosmic neutrinos manifest themselves directly in the process of forming the lightest nuclei 
in the early universe so that Big Bang Nucleosynthesis (BBN) yields compelling evidence 
for the reality of the "cosmic neutrino sea." 



4.4 Big Bang Nucleosynthesis 33 

All of the deuterium and most of the helium in the present-day universe can not have 
been produced in stars, so that it must have been produced shortly after the Big Bang. The 
process of primordial nucleosynthesis explains today's abundance of several light elements, 
notably 4 He, D, 3 He and 7 Li. Perhaps the most important quantity is the mass fraction Y p 
of primordial 4 He, which is observationally inferred to be 

Y p = ^ Mu ° ~ 0.24. (85) 
p M H + M He V ; 

The abundances of the other elements are very much smaller; for example ^D+SHe/^H — 
1(T 5 , n7 U /n u ~ 1(T 10 . 

In a thermal plasma, all nuclei should be present with their thermal equilibrium abun- 
dance given by nuclear statistical equilibrium (NSE). The abundance of nuclei with atomic 
number A, nuclear charge Z, binding energy Ba, and a statistical factor qa is || 

n A = 9 A A^2~ A [£y) 2 n^niT z expB A /T. (86) 

For example, Qa — 2 for 3 He which has total spin 1/2. 

At high temperatures, there are very few nuclei because the dissociated state is pre- 
ferred. The early universe is a "high-entropy environment" because there are about 10 10 
thermal photons per baryon. Therefore, the dissociated state is favored until very late when 
the temperature is significantly below the typical nuclear binding energies. The main idea 
behind BBN is that nuclei stayed dissociated until very late so that they appear in appre- 
ciable numbers only at a time when the reactions forming and dissociating them begin to 
freeze out of thermal equilibrium. Therefore, hardly any nuclei heavier than helium form 
at all — otherwise one might have expected that most of the matter appears in the form 
of iron, the most tightly bound nucleus. 

Nuclei are produced from protons and neutrons which, at high temperatures, are in 
thermal equilibrium by /9-processes of the form 

v e + V e+ + n i v e + n <-> + p, n <-> p + e~ + v e . (87) 

These processes "freeze out" when they become slow relative to the expansion rate of the 
universe. By dimensional analysis the reaction rates are of the form 

T oc G 2 F T 5 . (88) 

The other characteristic time scale is defined by the Hubble parameter H which has the 
form 

H oc — T 2 , (89) 
m P i 



34 4 Neutrinos in Cosmology 



where we write Newton's constant as the inverse of the Planck mass squared. Therefore 
the expansion of the universe was too fast for the rates when H > T, which leads to a 
freeze-out temperature of 

T F ~ 0.8 MeV. (90) 

At lower temperatures (later times) weak processes and thus neutrinos are no longer in 
equilibrium. 

Before the weak reactions became ineffective, the ratio of protons and neutrons was 
given thermal equilibrium as 



^ = exp ( -#) , (91) 



n p \ T 

where Q = m n — m p = 1.293 MeV. When T < Tp (after about 1 s), the ratio changes only 
by neutron free decays with a neutron lifetime of r n ~ 887 s, i.e. it decreases as exp (—t/r n ). 

At Tp when the weak interactions froze out, the neutron/proton ratio was about 1/6. 
After the temperature fell below about 0.1 MeV at a cosmic age of about 3 min, BBN 
began in earnest. Some important reactions involving the first few steps are 

n + p «-> D + 7 (E B = 2.2 MeV), 
D + p^ 3 He + n, D + n <-> 3 H + 7, D + D <-> 3 He + n, (92) 

3 He + n <-> 4 He + 7, 3 H + p <-> 4 He + 7, D + D <-> 4 He + 7. 

Most of the reactions end in 4 He. Heavier elements do not efficiently form because they 
have too large Coloumb barriers. Moreover, there are bottlenecks in the reaction network 
because there are no stable nuclei with A = 5 or A = 8. 

At higher temperatures the deuteron produced in the first reaction was immediately 
dissociated by photons, which have an average energy of (E 7 ) = p 7 /n 7 ~ 2.70 T. These 
dissociating photons were available in great number since the baryon/photon ratio 77 = 
?7iolO~ 10 is very small. When r\ exp (Eb/T) ~ 1, deuterons were produced in sufficient 
numbers. This requirement gives a temperature of about 0.1 MeV and thus defines the 
beginning of BBN. 

When all remaining neutrons had disappeared in 4 He, we have Mn e = n n 1m^ and 
Mu = (n p — n n )mN, which leads to a BBN prediction for the helium mass fraction of 

Y p = (93) 

Tip T" Tl n 

Since BBN stopped about three minutes after the weak reactions Eq. ( |57D froze out, we 
get n n /rip ~ 1/7 and thus Y p ~ 0.24 in agreement with the observations. 

All of our estimates depend on the freeze-out temperature Tp. It would be higher if 
H were larger than implied by the standard value for g*. A larger Tp implies an increased 
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Figure 8: Typical BBN predictions (solid line) for^,, D, 3 He and 7 Li relative to H with 
2cr theoretical uncertainty. The vertical band represents a recent measurement of the 
deuterium abundance []3~^| . (Figure from |33|| .) 
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1/2 

n n /n p and thus an increase of Y p . Since H is proportional to gj we find that additional 
neutrino families would lead to a higher Y p . This argument was used to constrain the 
number of neutrino flavors before the famous LEP experiment at CERN. In [|33] one finds 
that for Y p = 0.246 ± 0.0014 the number of neutrino species is limited by N v < 3.20 at 
95% C.L. 

The theoretical dependence of the light element abundances on rj is shown in Fig. [8], 
assuming the standard number of neutrino families N v = 3. The BBN predictions explain 
abundances consistently over a range of 10 orders of magnitude. Therefore, at the present 
time it appears that BBN is a consistent theory. It certainly requires the presence of 
thermal neutrinos and thus can be taken as compelling evidence for the reality of the 
cosmic neutrino sea. 



4.5 Neutrinos as Dark Matter 

The experimental evidence for neutrino oscillations is now so compelling that the reality of 
non-vanishing neutrino masses must be taken as a serious working hypothesis. However, 
neutrinos are not a good candidate for the ubiquitous dark matter (DM) which dominates 
the dynamics of the universe. 

One of the most conspicuous dark matter problems is that of the flat galactic rotation 
curves, implying that the dynamical mass of the galaxy is dominated by some non-luminous 
component. Assuming that this mass density consist of neutrinos, their maximum phase- 
space density implied by the Pauli exclusion principle is 

Here, p max = m v v max is the maximum momentum of a neutrino bound to the galaxy with 
■^max the escape velocity of order 500 km s _1 . Since p v xm — n u m v we obtain a lower limit 
for m v . For typical spiral galaxies one finds m v > 20-30 eV, and much larger values for 
dwarf galaxies which are dominated by dark matter. This limit, known as the Tremaine- 



Gunn bound |34| , at least nominally excludes neutrinos as the main component of galactic 
dark matter. 

Even more severe problems arise from arguments about cosmic structure formation. 
The universe is thought to have started from an almost homogeneous early phase with 
low-amplitude primordial density fluctuations which must have been produced by some 
physical process; one favored scenario involves the generation of the primordial fluctuation 
spectrum during an early phase of exponential expansion (inflation). The density contrast 
of these primordial fluctuations increases by the action of gravity since any region with 
more mass than its surroundings will attract more mass, at the expense of lower-density 
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regions. This gravitational instability mechanism can beautifully account for the observed 
structure in the matter distribution. 

However, weakly interacting particles become collisionless early and thus can travel 
large distances undisturbed. In this way they will wash out the primordial density fields 
by their "free streaming" or "collisionless phase mixing" and thus remove the seeds for 
structure formation up to a certain scale. The smaller the mass, the later these particles 
will become nonrelativistic, the further they travel, and thus the larger the scales below 
which the primordial field of density fluctuations has been erased. If all scales up to those 
which later correspond to galaxies have been erased, the particles are called "hot dark 
matter" while otherwise they are referred to as "cold dark matter." The division line 
between the two notions corresponds to a particle mass in the few keV range. 

Since hot dark matter can not account for the observed structures within the standard 
model of structure formation, they are strongly disfavored, except perhaps as a subdom- 
inant component in a hot plus cold dark matter scenario. It is quite surprising that 
oscillation experiments have established neutrino masses, but that neutrinos nevertheless 
are no good dark matter candidates. 




m v [eV] 

Figure 9: Contribution of massive neutrinos to Q. 

It is widely believed that the dark matter consists of some new, heavy, neutrino-like 
particles generically referred to as WIMPs (weakly interacting massive particles). One may 
think that such a cold dark matter particle would vastly overdose the universe since the 
cosmological mass limit on neutrinos requires a mass less than a few ten eV. However, if 
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the weakly interacting particles are heavier than a few MeV, the weak-interaction freeze- 
out happens when the temperature has fallen below their mass. Therefore, their number 
density is decimated by annihilations, i.e. their thermal equilibrium density is suppressed 
by a Boltzmann factor e~ m l T until the annihilation rate becomes too slow to keep up 
with the cosmic expansion. Therefore, as a function of the assumed mass of a weakly 
interacting particle one obtains a cosmic density contribution schematically shown in Fig. [|. 
This "Lee-Weinberg curve" turns over for large neutrino masses so that there is a second 
solution at masses of a few or a few ten GeV where neutrinos could be the dark matter, 
and then of the cold variety. 

Of course, none of the sequential neutrinos can play this role because their masses are 
too small. Even if we ignore the oscillation data, the direct kinematical mass limits are 
so restrictive (about 18 MeV for u T , the worst case) that standard neutrinos can not play 
the role of cold dark matter. However, in popular extensions of the standard model the 
existence of the requisite particles is quite plausible. Notably the theory of supersymmetry 
provides an ideal candidate in the guise of the "neutralino," a neutral spin-1/2 Majorana 
fermion. From a cosmological perspective, these neutralinos are virtually identical with a 
heavy Majorana neutrino. Even if the standard neutrinos are not the dark matter of the 
universe, it still looks like weakly interacting particles make up the bulk of the mass in our 
universe. 
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5 Conclusions 

The question if neutrinos, the most elusive of all elementary particles, have nonvanishing 
masses has come very close to a conclusion over the past few years. For a long time it 
had been suspected that the solar neutrino problem is solved by neutrino oscillations and 
thus by neutrino masses. However, it took the spectacular up-down asymmetry of the 
atmospheric neutrino flux observed by the SuperKamiokande detector to convince most 
of the sceptics that neutrino oscillations are real. It is remarkable that natural neutrino 
sources, the Sun and the upper atmosphere, play a central role at putting together the 
pieces of the jigsaw puzzle of the leptonic CKM matrix. 

As we learn more about the neutrino masses and mixing parameters, there is a certain 
sense that neutrino masses are too small to be of great cosmological importance. To be 
sure, it is still possible that the global neutrino mass scale is much larger than their mass 
differences. In such a degenerate scheme, neutrinos could still play a role as a subdominant 
dark matter component in hot plus cold dark matter cosmologies. It is also still possible 
that sterile neutrinos exist. If this hypothesis were verified, neutrinos would be assured of 
an important role for cosmic structure formation and might modify the standard theory of 
Big Bang nucleosynthesis. 

One should take note that neutrino astrophysics is a field much broader than the ques- 
tion of the neutrino mixing matrix. Once the mixing parameters and masses have been 
settled, one has greater confidence in one's treatment of astrophysical phenomena were 
neutrino masses and oscillations are potentially important such as supernovae or gamma- 
ray bursts, or in the use of neutrinos as a new form of radiation to do astronomy with. 
Neutrino physics and astrophysics will keep us busy for some time to come! 
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A Useful Tables 
A.l Integrals 

The following integrals frequently appear in the context of calculations involving particle 
reactions in thermal media, where ( refers to the Riemann zeta function. 

Table 1: Thermal integrals. 
Maxwell-Boltzmann Fermi-Dirac Bose-Einstein 

oo oo oo 

— dx / dx / dx 

e x J e x + 1 J e x -1 

o oo 

n = 2 2 |Cs ^ 1-8031 2( 3 ^ 2.40411 

n = 3 6 ~ 5.6822 £ ~ 6.4939 



A. 2 Conversion of Units 

We always use natural units where % = c = k-Q = 1. In order to convert between different 
measures of length, time, mass, or energy one may use the following table. For example, 
1 cm -1 = 2.998 x 10 10 s -1 . The atomic mass unit is denoted by amu. 



Table 2: Conversion factors for natural units. 





s- 1 


cm 1 


K 


eV 


amu 


erg 


g 


s- 1 


1 


0.334xl0~ 10 


0.764xl0- n 


0.658xl0~ 15 


0.707xl0~ 24 


1.055xl0~ 27 


1.173xl0~ 48 


cm -1 


2.998xl0 10 


1 


0.229 


1.973xl0~ 5 


2.118xl0~ 14 


3.161xl0~ 17 


0.352 xlO" 37 


K 


1.310xlO n 


4.369 


1 


0.862xl0~ 4 


0.962 xl0~ 13 


1.381X10" 16 


1.537xl0~ 37 


eV 


1.519xl0 15 


0.507xl0 5 


1.160xl0 4 


1 


1.074xl0~ 9 


1.602xl0~ 12 


1.783 xlO" 33 


amu 


1.415xl0 24 


0.472xl0 14 


1.081 xlO 13 


0.931 xlO 9 


1 


1.492xl0" 3 


1.661 xl0~ 24 


erg 


0.948xl0 27 


0.316xl0 17 


0.724xl0 16 


0.624xl0 12 


0.670xl0 3 


1 


1.113xl0- 21 


g 


0.852 xlO 48 


2.843 xlO 37 


0.651 xlO 37 


0.561 xlO 33 


0.602 xlO 24 


0.899 xlO 21 


1 
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